DESY 12-007 
IFUM-986-FT 



(N 



Threshold resummation in SCET vs. perturbative QCD: 

an analytic comparison 

Marco Bonvini, 1 ' 2 * Stefano Forte, 3 Margherita Ghezzi 3 ' 4 ^! and Giovanni Ridolfi 1 



1 Dipartimento di Fisica, Universita di Genova and IN FN, Sezione di Genova, 
Via Dodecaneso 33, 1-16146 Genova, Italy 
■ 2 Deutsches Elektronen-Synchroton, DESY, 

. Notkestrafie 85, D-22603 Hamburg, Germany 

3 Dipartimento di Fisica, Universita di Milano and INFN, Sezione di Milano, 
5"! Via Celoria 16, 1-20133 Milano, Italy 

4 Dipartimento di Fisica, Sapienza Universita di Roma and INFN, Sezione di Roma, 
^ ; Piazzale Aldo Moro 2, 1-00185 Roma, Italy 

CM 

£^ Abstract: 
Si- 

We compare threshold resummation in QCD, as performed using soft-collinear effec- 
. tive theory (SCET) in the Becher-Neubert approach, to the standard perturbative QCD 

formalism based on factorization and resummation of Mellin moments of partonic cross- 
sections. We consider various forms of the SCET result, which correspond to different 
choices of the soft scale fi s that characterizes this approach. We derive a master formula 
^sO I that relates the SCET resummation to the QCD result for any choice of fx s . We then 

H ' use it first, to show that if SCET resummation is performed in iV-Meflin moment space 

by suitable choice of fj, s it is equivalent to the standard perturbative approach. Next, we 
show that if SCET resummation is performed by choosing for /i s a partonic momentum 
variable, the perturbative result for partonic resummed cross-sections is again reproduced, 
but like its standard perturbative counterpart it is beset by divergent behaviour at the 
endpoint. Finally, using the master formula we show that when [i s is chosen as a hadronic 
momentum variable the SCET and standard approach are related through a multiplicative 
(convolutive) factor, which contains the dependence on the Landau pole and associated 
divergence. This factor depends on the luminosity in a non-universal way; it lowers by 
one power of log the accuracy of the resummed result, but it is otherwise subleading if 
one assumes the luminosity not to contain logarithmically enhanced terms. Therefore, the 
SCET approach can be turned into a prescription to remove the Landau pole from the 
perturbative result, but the price to pay for this is the reduction by one logarithmic power 
of the accuracy at each order and the need to make assumptions on the parton luminosity. 
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1 Threshold resummation and the Landau pole 



The interest in the resummation of logarithmically enhanced contributions due to soft 
gluon radiation in perturbative QCD (threshold resummation, henceforth) has been re- 
cently revived due to its relevance for many LHC processes, such as Higgs [1] or top [2] pro- 
duction. Threshold resummation was originally performed (to next-to-leading logarithmic 
accuracy) by factorizing the hadronic cross-section in Mellin space in terms of a luminosity 
and a partonic cross-section, and then exponentiating logarithmically enhanced corrections 
to the latter to all orders through eikonal [3] or factorization [I] techniques. Subsequent 
derivations and generalizations to all logarithmic orders were obtained, among others, from 
a suitable two-scale generalized factorization theorem [5] and through renormalization- 
group improvement of the kinematics of the gluon radiation phase-space [6], with an 
additional hypothesis of factorization of virtual corrections. 

In all these approaches, resummation is performed after Mellin transformation of the 
hadronic cross-section, which factorizes it into the product of a parton luminosity and 
a partonic cross-section. More importantly, in Mellin space the partonic cross-section in 
the soft limit can be obtained by exponentiating single-particle emission cross-sections, 
thanks to the fact that in Mellin space the n-particle longitudinal phase space factorizes. 
The large logs which are resummed are then logarithms of N, the variable which is Mellin 
conjugate to r (a dimensionless ratio which equals one at threshold), rather than the 
original ln(l — r). 

More recently, factorization and exponentiation were directly performed at the level 
of Feynman diagrams, without the need for a Mellin transformation, using path-integral 
methods to separate off soft gluon modes [TIE] - I n the latter approach the standard re- 
summed results are readily recovered, but the way the terms which dominate in the eikonal 
limit emerge order-by-order in perturbation theory (and the next-to-eikonal corrections to 
them) is particularly transparent. Indeed, an important use of resummed results is to pro- 
vide predictions for higher order terms which can even be used to construct approximate 
expressions for unknown fixed-order corrections (see e.g. Ref. [S]). 

However, regardless of how resummation is proven, resummed expressions for partonic 
cross-sections with a fixed logarithmic accuracy in momentum space (i.e. next fc -to-leading 
ln(l — z), where z is a partonic scaling variable) turn out to be ill-defined: they lead 
to divergent hadronic cross-sections upon convolution with a parton luminosity [6jll0|. 
This behaviour is already present at the fixed-coupling level |10| . and it persists when 
the coupling runs [6]. It can be traced [10] to the fact that the truncation of resummed 
results to any finite logarithmic accuracy in momentum space induces terms which vio- 
late longitudinal momentum conservation, thereby leading to factorial divergence of the 
perturbative expansion: the result is well-defined provided only the truncation to finite 
logarithmic accuracy is performed in Mellin space (i.e., next fc -to-leading lniV, rather than 
next fc -to-leading ln(l — z)), and the Mellin transform is inverted exactly to power accuracy, 
i.e. retaining terms to all logarithmic orders in 1 — z and only neglecting terms which are 
down by powers of 1 — z [6J. As a consequence, perturbative QCD resummation, even if 
derived using a momentum-space argument, must be performed in Mellin space (to finite 
logarithmic accuracy) if it is to respect momentum conservation, and to lead to finite 
physical (hadronic) cross-sections. 

At the running-coupling level, however, a new difficulty arises: namely, it turns out 
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that the next -to- leading In ./V series of contributions to the partonic cross-section at any 
finite logarithmic accuracy, viewed as a series in the strong coupling a s , corresponds, 
upon inverse Mellin transformation, to a divergent series of contributions to the partonic 
cross-section. This divergence can be traced to the Landau pole in the strong coupling: 
as long known [llj . resummed results correspond to effectively replacing the hard scale 
M 2 at which the strong coupling is evaluated with a scale M 2 (l — z) a related to the soft- 
gluon radiation process (with a a process-dependent exponent, e.g. a = 2 for Drell-Yan). 
Because the hadronic observable is found by convoluting the partonic cross-section with a 
luminosity , the integration over parton momenta always intercepts the region z — > 1 where 
the strong coupling blows up, and this manifests itself as a divergence of the expansion in 
powers of a s (M 2 ). This divergence, which is of non-perturbative origin, can be removed 
by addition of subleading terms: within the commonly used "minimal prescription" of 
Ref. |10] this is done by choosing a particular integration path to perform the Mellin 
inversion integral, which corresponds to adding a term which is more suppressed than any 
power of 1/M 2 , while with the more recent "Borel prescription" |12tll3| this is done by 
adding a higher twist term to make the divergent series Borel summable. 

An alternative approach to resummation can be based on the soft-collinear effective 
field theory (SCET) [T3HTT] . which provides [18] an alternative derivation of QCD fac- 
torization: threshold resummation based on SCET was performed in Refs. [19H23] . This 
approach provides a powerful alternative way of determining resummed results for hadronic 
observables, which can then be used for phenomenology through the standard Mellin-space 
formalism of Ref. [10]. However, it was pointed out in Ref. [23] that, thanks to the fact 
that the effective theory deals with the hadronic degrees of freedom, in a SCET approach 
resummed expressions can be directly derived in terms of the hadronic kinematic variable, 
i.e., in practice, SCET allows one to perform the resummation of ln(l — r), where r is a 
measurable dimensionless kinematic ratio. The advantage is that the divergences related 
to the need to integrate over the parton kinematics are no longer present: hence, in par- 
ticular, the difficulties related to the Landau pole of the strong coupling disappear. The 
approach of Ref. [23] has been subsequently developed for phenomenology, and applied 
to various physical processes, such as deep-inelastic scattering [24] . Drell-Yan [25J and 
Higgs [26] production. 

Henceforth, for brevity, we will refer to the approach of Ref. |23j as SCET approach, 
and that of Refs. [3H6] as QCD approach. It should be observed, however, that whereas 
the QCD implementation of resummation is unique to any finite perturbative order, the 
aforementioned |19H23| alternative implementations of threshold resummation in SCET, 
to the best of our understanding, lead to results which differ (possibly by subleading terms) 
even when truncated to finite perturbative order. Here we will concentrate on the SCET 
approach of Ref. [23], which has been widely used in particular for phenomenological 
applications. 

However, results obtained in the approach of Ref. |23] are not easily compared to 
those obtained using the standard approach of Refs. [3H6], because the direct connection 
to factorization and resummation at the level of partonic cross-sections is lost. Indeed, 
as mentioned, the presence of the Landau pole implies that the expansion in powers of 
a s (M 2 ) of the resummed partonic cross-section diverges. Hence, if the resummed SCET 
result of Ref. [23] is free of divergences, its expansion to fixed order must necessarily differ 
from that of the standard Mellin-space resummation. 
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This difference has never been determined so far: its computation is the goal of this 
paper. Clearly, its knowledge is crucial in order to determine the theoretical and phe- 
nomenological viability of the SCET resummation of Ref. |23j . Some phenomenological 
comparisons of resummed predictions for relevant physical processes obtained using SCET 
to standard perturbative results have been performed in Refs. |24t426j . Differences are 
found to be reasonably small: however, this does not shed light on their analytic form. 
But knowledge of this analytic form is necessary if we wish to know, first, whether up to 
the stated accuracy the SCET and QCD approaches are equivalent, and second, even if 
they are, what is the kind of subleading suppression of the terms introduced in the SCET 
approach to tame the perturbative divergence, i.e. what is the accuracy of the SCET 
approach (be it power or logarithmic). 

The answer to these questions is presented here in several steps. In Section [21 after 
summarizing the known form of resummed results both in the perturbative QCD and 
SCET approach, we recall that the definition of next fc -to-leading log accuracy in the SCET 
approach of Ref. [23] and in the standard perturbative QCD approach are different, and 
only agree at the leading log level. Beyond the leading log, SCET results are always less 
accurate by one power of log than the perturbative ones: so N fc LL in the perturbative 
case always include terms which only appear in the N fc+1 LL SCET result, and so forth. In 
order to proceed to a comparison, it is necessary to discuss the dependence of the SCET 
resummation on the soft scale: in Section [3] we summarize how SCET results in Mellin 
space, or in momentum space, at either the partonic or hadronic level can be obtained by 
different choices of soft scale. A comparison is then made possible through the derivation 
of a general relation between the SCET result and the standard result, by expressing the 
latter in terms of the convolution of the former with a function C r which depends on 
the soft scale. We establish this result at next-to-next-to-leading logarithmic order (but 
we conjecture it to hold to all orders): it provides a master formula which enables a full 
comparison of the QCD and SCET results, both from an analytic and a numerical point 
of view. 

As a preliminary step, this master formula can be used to prove the fact that if SCET 
resummation is performed in Mellin space it is completely equivalent to the standard 
approach, and in particular it has the same logarithmic accuracy at each order. This 
result was established already in Refs. [241125] . but with the aforementioned lower log 
accuracy of the SCET results. This is done in Section where we also digress to show 
that if SCET resummation is performed in momentum space by choosing a partonic scaling 
variable z, it coincides with the perturbative result up to power suppressed corrections, 
but, like the perturbative result, it diverges at the partonic endpoint z = 1. We can 
then (in Section \5§ tackle the computation of the function C r which relates the SCET 
and perturbative resummation when the soft scale is chosen as a measurable hadronic 
scale. In this case, the SCET result is free of Landau pole, and thus the divergence is 
entirely contained in the C r function. This function depends on the PDF luminosity in a 
non-universal way, and thus whether or not it is subleading depends on the form of the 
PDF. In particular, if one assumes that the luminosity does not contain logarithmically 
enhanced terms, then we can show that this function is always logarithmically subleading, 
provided only the less accurate SCET definition of logarithmic accuracy is used. However, 
any logarithmically enhanced contribution to the parton luminosity C{x) proportional to 
ln fc (l — x), with k > 1, will lead to contributions to C r which are of the same order as 
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those induced by perturbative resummation. 

Therefore, we conclude that it is only for a particular class of luminosities that SCET 
with a hadronic choice of soft scale reproduces the perturbative result, and can thus 
be considered to be equivalent to the standard approach and to provide an alternative 
prescription to remove the divergence of the perturbative expansion. Even when this 
is the case, the momentum-space SCET resummation prescription of Ref. |23j requires 
lowering by one order (one power of log) the accuracy of the resummed result at each 
logarithmic order. Furthermore, in the SCET prescription, terms which are introduced in 
order to remove the perturbative divergence are only logarithmically subleading, rather 
than being power suppressed (as in the Borel prescription) or exponentially suppressed 
(as in the minimal prescription), along with power suppressed terms. Finally, subleading 
terms which are induced by SCET resummation are suppressed by powers or logs of the 
hadronic scale: this feature of SCET resummation may also be a limitation, because the 
partonic and hadronic scales, though related, do not coincide, and in fact it may well be 
that the former is close to threshold while the latter is not |27j . 



2 Threshold resummation at fixed logarithmic accuracy 

For definiteness, we will concentrate on the production of Drell-Yan pairs at hadron col- 
liders. This choice does not entail loss of generality, and the extension to other processes 
is straightforward. We will consider in particular the invariant mass distribution -SSz , 
with M the invariant mass of the pair. We define the hadronic scaling variable 

M 2 , , 

t = — (2.1) 

where s is the hadronic center-of-mass energy squared, so the threshold limit is r — > 1. 
Perturbative QCD factorization takes the form 

a(r,M 2 ) = £^C(z,M 2 )£g), (2.2) 

where C is the parton luminosity, and a(r,M 2 ) is a dimensionless cross-section 

defined by requiring that at the Born level (i.e. at order a®) C(z, M 2 ) = 5(1 — z). Note that 
Eq. (|2.2p is a schematic expression: in general, a sum over different parton subprocesses 
must be included. In the sequel, without significant loss of generality, we shall always 
choose the renormalization and factorization scales equal to each other and to the physical 
hard scale fi 2 F = ii 2 R = M 2 . 

2.1 Perturbative QCD: resummation in TV space 

As discussed in Section [H standard QCD resummation is more conveniently performed by 
taking a Mellin transform 

a(N, M 2 )= f dr t n - V(r, M 2 ) = C(N, M 2 )C{N) (2.4) 
Jo 
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which factorizes both the convolution Eq. (|2.2p and the gluon radiation phase space. In 
Eq. by slight abuse of notation we denote with C(N,M 2 ) and C(N) the Mellin 

transforms of C(z,M 2 ) and C(z) respectively. 

The iV-space resummed coefficient function has the form 

- / M 2 \ 

C QCD (N,M 2 ) = g (a s ) expS ( M 2 , —r ] (2.5) 



where 



M 2 \ 



M 2 , -z-TT | = / dzz N ~ 



N 2 J 



N 2 J 

1 rM^i-z? d[l 2 



1 



f Z d J^2A (a s (^ 2 )) + D (a s ([l - z] 2 M 2 )) 
Jm 2 H 



(2.6) 

The functions go(a s ), A(a s ) and D(a s ) are given as power series in a s , with go(0) = 1 
and A(Q) = D(0) = 0; A(a s ) is order by order the coefficient of the soft singularity in the 
Altarelli-Parisi splitting function for the relevant partonic subprocess, while the functions 
D(a s ) and go(a s ) are process-dependent. Specifically, in the case of Drell-Yan production 
initiated by quark-antiquark collisions, the relevant Altarelli-Parisi splitting function is 

P«(ol.,x) = [1 + 0(1 - *)] • (2.7) 

As a result, the resummed coefficient function takes the form (using the notation of 
Ref. [ID]) 

C QCD (Ar,M 2 ) = 5o (as)exp5(aL,a), (2.8) 

S(aL, a) = —gx(aL) + g2(aL) + ags(aL) + a 2 gAaL) + . . . , (2-9) 
a 

a = 2a s (M 2 )f3 , L = ln^, (2.10) 
where /3q is the first coefficient of the QCD f3 function, defined as 

S^ = -^) + 0(al>, h= UCA 1 - 2n ' (2.11) 

and the functions g^ are of order g\{aV) = 0(a 2 ) and gi(a>L) = 0(a s ) for i > 1, and are 
straightforwardly obtained performing the integrals in Eq. (j2.6fl . and thus each determined 
by a finite number of coefficients in the expansion of the functions A and D. Note that the 
functions go and S do not coincide with <7o and <S of Eq. (|2.5|) . because, by definition, S 
unlike S does not contain terms which are not logarithmically enhanced, while go includes 
non- logarithmic contributions both from go itself, and from the integral Eq. (|2.6|) . 

The standard perturbative QCD resummation predicts correctly all contributions to 
In C QCD (A / ", M 2 ) up to a given order: in other words, if contributions up to g n are included 
in S(aL,a) Eq. (|2.9p then In C QCD (N, M 2 ) is determined up to subleading corrections of 
order 0{a^ + ^ n 1 ^L fc ). This is standardly called N n_1 LL resummation. However, once 
the exponential is expanded out in order to obtain the coefficient function C QCD (A^, M 2 ), 
at each order in a s , only a restricted number of logarithmically enhanced terms is pre- 
dicted correctly, and furthermore, inclusion of the prefactor go(a s ) Eq. (|2,8p (which is not 
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logarithmically enhanced) is mandatory in order to improve the accuracy beyond NLL. 
In fact its inclusion already at the NLL level increases the number of contributions to 
the coefficient function which are predicted correctly. In Tab. [1] we summarize the order 
up to which the expansion of the functions Eq. (|2.8I2,9|) should be included to achieve 
a given logarithmic accuracy and, in the last column, the order of the contributions to 
the resummed coefficient function C QCD (iV, M 2 ) which, as a consequence, are predicted 
correctly. 

log approx. gi up to go up to order accuracy: a"L fc 

LL i = l (a s )° k = 2n 

NLL % = 2 (as) 1 2n - 2 < k < 2n 

NNLL i = 3 (a s ) 2 2n - 4 < k < 2n 



Table 1: Orders of logarithmic approximations and accuracy of the predicted logarithms 
in perturbative QCD. 



2.2 The SCET approach 

Resummation in SCET in the approach of Ref. [23], which henceforth we will refer to 
as SCET resummation for short, is directly given in the physical space of momentum 
fractions. The relevant expression for Drell-Yan pair production has been computed in 
Ref. [25], and it is given by 

C SCET (z, M 2 , nl) = H(M 2 )U(M 2 ,fi 2 )S(z, M 2 , /x 2 ) (2.12) 

where H(M 2 ), the so-called hard function, has an expansion in powers of a s computed at 
the hard scale M 2 ; 

/ jvf2 d \ 1 ( 1 - z \ 2r> e~ 2 ^ v 



where 



V= ^-Tcusp (a s (/x 2 )) ; T cusp (a s ) = A(a s ) (2.14) 
Jm 2 I 1 

and §dy(L,h) has a perturbative expansion in powers of a s (// 2 ). Note that the function 
r cU sp( a s) coincides with the function A(a s ) of Eq. (|2.6p . Finally, 



U(M 2 ,ti 2 s ) = exp | - / 

J A 



dfi 2 



M 2 V 2 



„2 

r cU s P (a s (/i 2 )) In ^ ~ iw («s(At 2 )) 



(2.15) 



where 7^(a s ) has a power expansion in a s . The resummed expression as given in Ref. [25] 
actually depends on several energy scales, which here for simplicity are all taken to be equal 
to the hard scale M 2 . 

Two important formal aspects characterize the SCET resummed result. The first is 
that it depends on a "soft scale" fj, s , and in fact the logs which are being resummed in 
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SCET are In ^f-. Hence, different choices of soft scale lead to different forms of the SCET 
resummation, as we shall discuss in greater detail in the next Section. 

The second is related to the well-known fact that at the endpoint z = 1 the coeffi- 
cient function C SCET (z, M 2 , fj, 2 ) is a distribution, rather than an ordinary function. This 
distribution is usually expressed in terms of the so-called plus distribution • The 

distributional nature of the SCET result emerges in the following way. The convolution 
product of C SCET (,z, M 2 , jit 2 ) with any well-behaved test function of z is well defined as 
long as 7] is a fixed, positive number: the factor (1 — z) 2r] acts as a regulator of the soft 
singularity at z = 1. The result can then be analytically continued to negative values of r? 
(which is typically the case in DY-like processes) by means of the identity 

dz (1 - z) 2 ^ 1 f(z) = f dz (1 - zf^ 1 [f(z) - /(l)] + ±f(l). (2.16) 
Jo 2r? 

Eq. (|2.16p defines a distribution on a space of test functions f(z), regular in the range 
< z < 1, which is usually written 

(1 - z) 2 "- 1 = [(1 - z)^- l ] + + ±.5(1 - z). (2.17) 

It is important to note that 77 is of order q s : therefore, the term proportional to 5(1 — z) 
in Eq. (|2.17|) combines with the factor 1/T(27?) = 2r] + 0(rj 2 ) in Eq. (|2.12|) to form an 
order-Og contribution (with the correct kinematic structure). 

As in the perturbative case, a given logarithmic accuracy is obtained by including a 
finite number of terms in the perturbative expansion of the functions which determine 
the resummed result, namely r cusp , jw, H and sdy- The accuracy which, according to 
Ref. [25] . is obtained by including in the SCET expression Eq. ([2.12H2.15]) coefficients up 
to a given order, as well as the corresponding nomenclature, are summarized in Tab. [2j As 
in the case of Tab. [TJ the last column gives the order of the contributions to C SCET which 
are predicted exactly. As mentioned in Section [1] and as is apparent comparing Tab. [Q to 
Tab. El beyond LL the SCET results are always less accurate than the QCD results of the 
same name: the QCD NLL includes terms of order a™ ln k jj with k > 2n — 2, but the 
SCET NLL only includes terms with k > 2n — 1. This was already observed in Ref. |28j . 



RG-impr. PT 


log. approx. 


A cusp 


lw 


H, sdy 


accuracy: a" ln fc ^ 




LL 


1-loop 


tree-level 


tree-level 


k = 2n 


LO 


NLL 


2-loop 


1-loop 


tree-level 


2n - 1 < k < 2n 


NLO 


NNLL 


3-loop 


2-loop 


1-loop 


2n - 3 < k < 2n 


NNLO 


NNNLL 


4-loop 


3-loop 


2-loop 


2n - 5 < k < 2n 



Table 2: Different approximation schemes for the evaluation of the resummed cross- 
section formulae in the SCET approach. 



When comparing the two different definitions of logarithmic accuracy, Tab. Q] and 
Tab. [21 one should distinguish a purely terminological issue and an issue of substance. 
The terminological issue is how each given accuracy is called: this is clearly immaterial. 
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The issue of substance is whether at (say) NLL the SCET expression Eq. (j2.12M2.15]) may 
be upgraded to the higher accuracy of the NLL QCD expression Eq. (|2.5H2.6|) (without 
having to resort to the yet more accurate NNLL SCET expression), and likewise at all 
subsequent logarithmic orders. We will show that the answer to this question depends on 
the choice of soft scale /j s ^ 



3 Choice of the soft scale and SCET-QCD comparison 

In the standard perturbative QCD approach to soft resummation, the energy scale which 
characterizes soft gluon emission is of the order of M(l — z): when the observed final 
state carries away a fraction z of the available partonic energy, the energy available for 
unobserved radiation is M(l — z), which is much smaller than M if z is close to 1. The fact 
that the scale involved is partonic has phenomenological implications: because the partonic 
center-of-mass energy is always smaller than the hadronic one, threshold resummation may 
be relevant even for processes which are relatively far from hadronic threshold, provided 
the parton luminosity is peaked for small values of the momentum fraction [27]. This 
indeeds is known to happen for Higgs production in gluon fusion at the LHC. |30]I31| 

In SCET resummation, however, one resums logs of the large ratio M//x s of the hard 
scale M to the soft scale fj, s , and various choices for the soft scale fj, s are possible: in 
particular, the choice which has been advocated in Refs. |23H26| . and which removes the 
problem of the Landau pole, consists of choosing for fj, s a scale which characterizes the 
(hadronic) physical process. 

If n s is chosen as a function of the partonic scaling variable z, then the resummed 
SCET partonic cross- section C SCET (z, M 2 , fj 2 ) Eq. Q2H2J) can be directly compared to the 
momentum-space perturbative QCD expression, which may be obtained by determining 
the inverse Mellin transform Cq CD (z, M 2 ) of the resummed iV-space expression Eq. (|2.5p . 
We will study this case in detail in the next Section. However, if fj, s is chosen as a function 
of the hadronic scaling variable r, the SCET and perturbative QCD resummed results must 
be compared at the level of physical cross-sections <7q CD (t, M 2 ) and cr SCET (r, M 2 ), which 
are respectively obtained substituting C QCD (z, M 2 ) or C SCET (z, M 2 , /i 2 ) in the factorized 
expression Eq. (|2.2[) , with some particular choice of soft scale fx s . 

It is important to understand that these different choices of soft scale lead to resummed 
predictions with different analytic structure. To see this, note that if the soft scale only 
depends on the parton momentum fraction z, then Eq. (|2.2p is a convolution, in the sense 
that upon Mellin transformation it factorizes according to Eq. (|2.4|) . This factorization is 
of course a necessary and sufficient condition for parton radiation to respect longitudinal 
momentum conservation. But if in Eq. (|2.2p the coefficient function depends on r through 
the soft scale, then the convolution structure is destroyed. This means that with this 
particular choice of soft scale, upon Mellin transformation the cross-section no longer 
factorizes, thereby violating longitudinal momentum conservation. This also violates the 
QCD factorization theorem, because the short-distance partonic cross section depends on 
long-distance physics through the hadronic variable r. The possibility of making this 

In other contexts, such as for example the resummation of jet veto logs [29], SCET results which 
correspond either of two different accuracies, respectively akin to Tab. [5] or Tab. [T] may be achieved by 
suitable choices of terms to be included in the resummed expression. 
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choice of soft scale, and indeed the very possibility of making alternative choices of soft 
scale, some of which preserve factorization and some of which do not, appears puzzling 
in a standard perturbative QCD approach. We will not attempt to address the issue of 
principle of understanding this apparent structural discrepancy between SCET and QCD 
results. Rather, we will take the SCET and QCD expressions at face value: our aim will 
be to determine how they are related to each other. 

We will now derive a master formula which relates the SCET resummed expression 
for generic choice of the soft scale to the standard perturbative QCD expression. For def- 
initeness, we specialize to the next-to-next-to-leading log case, but all relevant structures 
are already present at this order so generalization to higher logarithmic orders is straight- 
forward. First, we give the explicit expression of the QCD result Eq. (|2.5p to this order. 
Then, we give the SCET expression Eq. (|2.12M2.15|) to the same order, and we perform 
its (exact) Mellin transform in order to allow for a comparison with the QCD expression, 
which is given in N space. Finally, by comparing the two expressions we derive a master 
formula which relates them, as a function of the soft scale fj, s , through a suitable factor 
(in Mellin space) or a convolutive function (in momentum space). 



3.1 Perturbative QCD resummation to NNLL 



The NNLL resummed expression in perturbative QCD is given by Eq. (|2.5p with [9ll32| 
(see also Ref. [27]) 



An 



A 



. - ^2 2 i J j, . , n , 



A 



3„,3 



4 
AC F 



16 



64 



7T 



D(a s ) = Dia s + D 2 a 2 s + 0(a 3 s ), 



D 1 = 0, D 2 



C F 
16vr 2 



C A 



1616 88 2 



27 



+ ¥ ^ + 56C 3 ) + ( — -^)n f 



224 16 2 



27 
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(3.1) 

(3.2) 
(3.3) 
(3.4) 



We can perform the z integral in Eq. ()2.5p using Eq. ()A.4j) : 



'QCD 



M 



M 2 

N 2 



M 2 



/i 



A (a s (u 2 )) ( In ^t- - In 



+ 



7T 



2 d 2 [ M2 ' N2 d^ 2 



12 dL 2 

2 



AP 



M 2 



C F 7T 

H ^— «s 



A(«s(/, 2 ))(lni,-lnA_ 



M 2 
N 2 



N 2 

2 N 1 



+ -D(a s (p 2 )) 



(3.5) 



where (as per Eq. (|A.5|) ) N = Ne 1 . We have neglected subleading (N 3 LL) terms (including 
the replacement N — > N in the argument of a s in the last term) and we have brought all 
integrals to a common form using 



dz 



1 - z 



M 2 (l-zf ,2 

%A{a s {» 2 )) 



M 2 



M 2 /N 2 du 2 f u 2 dfM 



M 2 



v Jm 2 A* 



A(a s (/i 2 )) 
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m 2 /n 2 dfj2 

M 2 A* 



\A{a a {^ 2 )) (in ^ " ln ^2 



1 - JZ 



£(a s (M 2 (l-z) 2 ) 



1 /• M2 / 7v2 d / u 2 



M 2 



-^D(a s ^ 2 )). 



(3.6) 
(3.7) 



In order to ease the subsequent comparison to the SCET result, we separate off the 
non-logarithmic constant from the last term in Eq. 



'qcd 



N 2 



M 2 /N 2 ^2 

M 2 V 2 
C F TT 



- l ("V)) (lni^-ln^) +D 2 a 2 s (v 2 ) 



+ 



a s (M' 



where 



- D 2 Cftt a Cf 



We can thus write 
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M 2 \ 



C QCU (N,M 2 ) = g (a s (M 2 ))e W S QCU ( M\j^\ 



where 



g (a s ) = 1 + g 01 a s + 0(a 2 ); 



M 2 



M 2 /N 2 



M 2 



M 2 



A (a s ( M 2 )) ( In - In ^ ) + D 2 a 2 a (ji 



N 2 



M 2 



,2/2\ 



(3.8) 
(3.9) 

(3.10) 

(3.11) 
(3.12) 



Note that go and S cannot be identified with go and 5 in Eq. (|2.5|) . because the integral 
in Eq. (|3.12p does contain some terms which are not logarithmically enhanced: 



A l f M2 ^ 2 dfi 2 



4 Jm 2 
so that 



I 1 



^a s (/i 2 ) ( ln-^ - In ^ ) = — 7 "as(M 2 )+logarithms+0(a 2 ), (3.13) 



N 2 



M 2 



2C F 



7T 



2C f o 

9oi = 9oi 7 ■ 



7T 



(3.14) 



However, the form Eq. (|3. 12|) of the exponent in the QCD result is especially suited for 
comparison to the SCET result, as we now show. 



3.2 SCET resummation to NNLL 

We turn to the SCET expression, which is given by Eq. (|2.12p with, to NNLL 



7w(u s ) = T w 



16vr 2 + 1 sj ' 



(2) 

Tw 



„ „ ( 808 llvr 2 \ ^ _ /224 4tt 2 \ 

C F C A (-- + — + 28 C 3 ) + C F T F n f \ — - - j . 



(3.15) 
(3.16) 
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In order to compare it to the perturbative QCD result, we perform a Mellin transform 
with respect to z. This is easy to do, because the z dependence is all contained in the soft 

whose Mellin transform is 



function S(z,M 2 ,nl 
M [S{z,M 2 ,n 2 s 



SDY 



, M 2 d 
Hi dr) 



-2777 



r(JV - tj)T(2tj) e 
T(N + V ) Tpq) 



Cp , o\ 



/#V 2 

It follows that the Mellin transform of the coefficient function Eq. (|2. 12jl is 

M 2 "l 



C SCET (iV,M 2 ,^) = tf(M 2 



1 + "^-<Wj + ^«s(ft, 



ln^ 



x cxp 



M 2 V 2 



cusp 



Of. 



.0* 2 )) In 



1 



In 



A* 

M 2 



+ 



(2) 

16vr 2 sVP 7 



+ 



(3.18) 



It is very important to observe that the Mellin transform has been computed at fixed 
fjL s . This means that firstly, Eq. (|3.18p is not the Mellin transform of the SCET expression 
when fi s depends on z (which we will discuss in the next Section): in that case the Mellin 
transform would also act on the z dependence through fj, s . And second, that if fj, s depends 
on r the cross-section o"scet(^~? M^ 2 ^) computed using Eq. (|2.2p does not factorize into the 
product of C SCET (N, M 2 , /j, 2 ) Eq. f|3. 18[) times a parton luminosity £(N) upon Mellin 
transformation: the Mellin integral over r would also act on the r dependence through fj, s 
which, as already noted, does not have the form of a convolution integral. 

Equation ()3.18j) can be brought in a form which is suitable for comparison to the QCD 
expression by separating off the constant as in Eq. (|3.8p , thus leading to 



C SCET (N,M 2 ,n 2 s ) = H(M 2 )E(N,M 2 ,n 2 s )expS SCET (M 2 ,^ 



with 



H{M 2 ) = H(M 2 )exp 
S SCET (M 2 ,fi 2 ) 



C F TT 

12 



a s (M 2 ) 



2~ 



.(2) 

Tw 



(2) 
TW 
16vr 2 



E(N,M 2 ,/i 2 



1 + 



M 2 A* 

C F TT 
12 

C F 



r cU s P («s (^ )) 



l + a s (M 2 ) 
1 



hi 



iV 2 



1% 



C F 
16vr 2 



2vr 



a s (/j, 2 ) ( In 



1 

N 2 



In 



Ms 
M 2 



(3.19) 

(3.20) 
(3.21) 
(3.22) 
(3.23) 



3.3 The master formula 

The QCD expression Eqs. (|3.1Uti3.12p and the SCET expression Eqs. (I3.19ti3.23p are easily 
related, by noting that, because r cusp (a s ) = A(a s ) and D2 = 7^ , the integrands in 
Eqs. (pmj) and (pl2Tp coincide, so 

5 S cet(M 2 , Q 2 ) = 5 QCD (M 2 , Q 2 ) = S(M 2 , Q 2 ). (3.24) 
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It follows that, splitting the integral as fj£ //P ^ = f$ 2 ^ + f*f //p we get 

C QCD (N,M 2 ) = C r (N,M 2 ^ 2 )C SCET (N,M 2 ^ 2 ) (3.25) 

where 



C r (N, M 



2 2 , 5o(« s (M 2 )) ex P 5 (/^^ 



(3.26) 



H(M 2 ) E{N,M 2 ,fj, 2 ) ■ 

The non-logarithmic terms in fact cancel to the accuracy of our computation. Indeed, 
by substituting the value [32] 



501 



C F 



IT 



(4C 2 - 4 + 2 7 2 ) 



in Eq. (|3.14p . and the value 



in Eq. ^20\i we get 



(3.27) 



(3.28) 



(3.29) 



ff(M 2 ) 

so deviations from unity are of the same order as the first contribution which, at NNLL 
accuracy, is not included in H(M 2 ) (according to Tab. [2]). The expression of C r can be 
further simplified by including the function E(N, M 2 , /j, 2 ) Eq. (|3.23p in the function S: 
indeed 



E(N, M 2 ,/i 2 ) = exp 



^a s (fi 2 s ) (in X- In 4> V 
8 s ^ s '\ jv 2 M 2 / 



+ 0(a 2 



exp 



? M 2 



/I 



-a s (fx 2 ) (in -i 



In 



iV 2 M 2 



—(3{a s (fi 2 )) f In - In At 



2 \ 2 



(3.30) 



Using Eq. (|3.30p in the definition of C r (N, M , /i 2 ) we obtain our final expression 



M 2 

C r (JV,M^) = expS r {/z*,-^ 



(3.31) 



with 



S r [i 



M 2 
N 2 



M 2 /N 2 ^2 

S M 2 



2 Aia s (/^ 2 ) 
^4(<W )) - - 



1 



+ ^i Uu ( / r))(ln^-ln-^ f ) + D,<K( fl 



4 

2 \ 2 



In — In -At 

N 2 M 2 



,2/,,2\ 



(3.32) 



Equation (|3.25|) together with the explicit expression Eqs. (|3.31H3.32"j) of the function C r 
provides the master formula which relates SCET and perturbative QCD resummation. It 
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is the main result of this paper. We note that no term of order a s appears in the integrand 
of Eq. (|3.32|) : indeed, the inclusion of the term E(N, M 2 , /i 2 ) has the effect of removing the 
term proportional to A\a s (see Eq. (|3.ip ). The remaining contributions to the integrand 
in Eq. (pT3"2"]) start at 0(a 2 ). 

It is important to observe that while Cq CD (N, M 2 ) does not admit a Mellin inverse, 
because it has a cut in the complex N plane starting at the value Nl at which the 
strong coupling blows up, C SCET (iV, M 2 , fi 2 ) does admit a Mellin inverse as long as fi s 
is kept fixed, because the argument of the strong coupling in the SCET expression does 
not depend on N. This means that if Eq. (|3.25|) is expanded in powers of a s (M 2 ), and 
then the expansion is Mellin-inverted term by term, the expansion of the left-hand side is 
divergent, while on the right-side the Mellin inverse of the expansion of C SCET (iV, M 2 , ji 2 ) 
converges to C SCET (z,M 2 ,n 2 ) Eq. ([2712]) . Therefore, the diver gence has been isolated in 
the Mellin inverse of the expansion of the function C r (N, M 2 , fi 2 ) Eq. (|3.26|) . 

If the perturbative expansion of both sides of Eq. (|3.25p in powers of a s (M 2 ) is trun- 
cated to any finite order, then the Mellin inverse of both sides exists, and one gets the 
momentum-space relation 

C QCD (z, M 2 ) = f ^-C r p, M 2 ,^ 2 ) C SCET (y, M 2 ,fi 2 ) , (3.33) 

where C SCET (z, M 2 , fi 2 ) is given by Eq. <|2. 12 j> (expanded out to the given order), while 
both C QCD (z, M 2 ) and C r (z, M 2 , // 2 ) should be understood as the truncation to the given 
order of the Mellin inverse of the expansion of the corresponding iV-space quantities. 
Equation (|3.33p is then the momentum-space version of the master formula. 

The master formula Eqs. (|3.25I3.33|) has been established at next-to-next-to- leading 
logarithmic order, defined according to Tab. [2j Note, however, that the accuracy is up- 
graded to the higher one of Tab. [I] if non-logarithmic terms cancel to 0(a 3 ), i.e. if the 
function H(M 2 ) in the SCET coefficient function Eq. (|3.19p is replaced by a function 
H(M 2 ) such that 

ffo(q s (M 2 )) = 3 
H{M 2 ) y sJ y ' 

Of course, this can always be achieved by letting H(M 2 ) = H(M 2 ) + H 2 a 2 (M 2 ) and 
suitably choosing the value of H2, while including the 0(a 2 ) to go(a 3 (M 2 )), as per Tab.[TJ 
(Whether H(M 2 ) coincides with the 0{a 2 ) expression of H(M 2 ) as obtained using SCET 
is an issue that we will not address here). We conclude that the master formula holds up 
to NNLL accuracy, defined as in Tab. [TJ It is easy to convince oneself that this argument 
should hold to all logarithmic orders. 

4 Perturbative QCD vs. SCET: partonic cross-sections 

The master formula Eqs. (j3.25H3.33]) shows how SCET resummation can be used to repro- 
duce standard results. Indeed, it immediately implies that if we fix the soft scale in terms 
of the Mellin-space variable, 

M , , 

M. = ^r, (4-1) 
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then C r (N,M 2 ,fi 2 s ) = 1, i.e. 



C QCO (N,M z ) = C SCET [N,M\ 



M 2 
N 2 



(4.2) 



so the standard QCD result is reproduced: with this choice, SCET resummation is per- 
formed at the level of Mellin-space partonic cross-sections. Notice that because with this 
choice the SCET and QCD expressions coincide, they also have the same accuracy. So 
with this choice the SCET results actually has the accuracy of Tab. [TJ rather than the 
lower accuracy of Tab.0 The equivalence of Mellin-space SCET resummation to the QCD 
expressions was already established in Ref. [33]; it was also pointed out in Ref. [25J, but 
with the lower accuracy of Tab. 

Alternatively, one may try to use SCET resummation for partonic cross-sections, but 
using the momentum-space SCET formula Eq. (|2. 12|) . with jj, s fixed as a momentum-space 
partonic scale, namely 

fi a = M(l-z). (4.3) 

This choice for instance was adopted recently in Ref. [M] to perform threshold resumma- 
tion for top production. This choice also provides another way of re-deriving the standard 
perturbative resummation from SCET. Indeed, it can be shown that, away from the end- 
point z = 1, all logarithmically enhanced terms ln A- ^ in the partonic cross-section are 
reproduced order by order with this choice. 

This is very easily seen at the leading-log, fixed-coupling level. Indeed, in this limit 
one has 

77 = ^iln(l-z), {1- z) 2r ' = exp[a s A 1 ln 2 (l- z)] , (4.4) 



so that 

C SCET (z,M 2 ,M 2 (l-z) 



cxp 



exp 



l-z) 2 " = exp [a s Ailn 2 (l-z)] , 



Aia 



M 2 



-^\n 2 (l-z) 



dp 1 £_ 

fj, 2 M 2 



{l-z) 2r i 1 
1-z T(2rj)' 



l-z r(2r?) 



(4.5) 



But to leading log order one may expand l/r(2r7) to first order in a s , so 



C SCET (z,M 2 ,M 2 (l- z) 2 ) = a s A 1 



ln(l - z) 
1 - z 



exp 



A\a s 



ln 2 (l-z) 



+ NLL; z^l. (4.6) 



On the other hand, the perturbative result in the same approximation is the inverse Mellin 
transform of 



C QCD (N,M 2 ) =exp 
i.e., using the results of Appendix IB"! 



a s A x 



In* I 

N 



NLL. 



C QCD (z,M 2 



l-z 



exp 



a s A x d 2 \ (1- z)Z 



2 aej no 



+ NLL; z^l. 



(4.7) 



(4.8) 



Expanding the exponential and keeping only leading log terms this is seen to coincide with 
Eq. flM]). 
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However, as pointed out in Ref. [10] and discussed in Section [IJ C QCD (z,M 2 ) (de- 
fined as the leading-log truncation of the inverse Mellin of Eq. (|4.7p ) is ill-defined at the 
endpoint z = 1: it behaves as a distribution which leads to a divergent integral upon 
convolution with any reasonably behaved luminosity, and, if expanded order by order in 
a s , it diverges factorially. The SCET expression, either in the form of Eq. (|4.5p or of 
Eq. (|4.6p . is also ill-defined as z — > 1. Indeed, because now rj depends on z (see Eq. (|4.4p ). 
it is no longer possible to use Eq. (|2.17p to regulate the behaviour of C SCET (z, M 2 , fi s ). 
Note that Eq. (|2.17p also had the effect of generating the required 0(a®) contribution to 
C SC et{z, M 2 , M 2 (I — z) 2 ) proportional to 5(1 — z). Furthermore, as z — > 1 the coefficient 
function Eq. (j4.5j) oscillates with a factorially-growing amplitude, because of the factor 
r(2T?) • ^he f ac t that ^ ne SCET resummed expression diverges at the partonic endpoint 
was already noticed in Ref. [31] ■ Because of these difficulties, we will not pursue further 
the choice Eq. (|4.3p of soft scale. 

5 Perturbative QCD vs. SCET in momentum space: hadronic 
cross-sections 

We now turn finally to the choice of soft scale which is recommended in Refs. |23H26| . 
specifically as a solution to the problem of the Landau pole, namely, a soft scale fixed in 
terms of the hadronic momentum scald 2 ! 



With this choice of soft scale, the SCET and perturbative QCD results can only be com- 
pared at the level of hadronic cross-sections 



Indeed, with the choice of soft scale Eq. (|5.ip the resummed SCET cross-section Eq. (|5.3p 
is no longer in the form of a convolution product, because the integrand depends on 
r explicitly in the lower integration bound and in the argument of C, but also implicitly 
through fi 2 . As a consequence, upon Mellin transformation with respect to r, <t scet (t, M 2 ), 
unlike the standard QCD result, does not factorize into a parton luminosity and a partonic 
cross-section. 

Therefore, the comparison must be carried out directly at the level of hadronic cross- 
sections Eqs. (|5.2H5.3p . using the momentum-space form Eq. (|3.33p of the master formula 
(always understood as a truncation to arbitrary but finite order in a s , as discussed in 
the end of Section [33]) . This is somewhat problematic, because the power counting of 
Tabs. [TH2] was defined at the level of coefficient functions and thus necessarily at the level 

In Refs. 23 -26 a slightly more general choice of soft scale is considered: namely, the soft scale Eq. (|5.1|l 
is generally rescaled by a function of r which does not vanish at r = 1, and is chosen in such a way that 
the finite-order perturbative expansion of sdy is reliable. Because this modification does not introduce 
any extra logarithmic enhancement, it does not affect our discussion, and we will not consider it. 



M» = M(l-r). 



(5.1) 




(5.2) 



(5.3) 
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of a partonic cross-section. Of course, it is possible to define a given logarithmic order 
at the level of SCET coefficient functions, then use this expression to compute the cross- 
section <r SCET (r, M 2 ) using Eq. (|5.3p . However, because this expression is not factorized, 
the question whether cr SCET (t,M 2 ) and o"qcd(t", M 2 ) agree at any given order can only 
be answered by comparing them directly, and counting logs of the hadronic scale 1 — r. 
The result will then inevitably depend on the choice of parton distributions. The only 
alternative is to simply conclude that the SCET result with this choice cannot be compared 
to the perturbative one, and cannot be endowed with a perturbative meaning [35J. 

We will perform this comparison by computing the difference between o" SCET (r, M 2 ) 
and 0" qcd (t, M 2 ) up to 0(a 2 (M 2 )) and using the master formula to relate results. We will 
then discuss the structure of the result to all orders. 

5.1 Fixed-order comparisons 

We start by computing the function C r (N, M 2 , /j, 2 ) Eq. (|3.3ip explicitly. Up to order a 2 
we find 



C(,Y. M-. f r) 1 ! n:i.U 2 ) ( -^A>ln 3 £ + ^ In 2 ± + 2D 2 In ± ) + 0(n*) (5. I ) 



where 

c = . (5.5) 

The corresponding momentum-space expression is readily obtained by performing the 
inverse Mellin transform of Eq. (|5.4[) with the help of Eq. ()A.9p : 

r2 ,.2\ 



C r (z,M 2 ,vi) = S(l-z) 



5=0 

(5.6) 



where the function 



K(z,0 = m^- 1 - z (5.7) 

plays the role of a generating function. 

The difference between the resummed physical cross-sections in the QCD and SCET 
formalisms is now found substituting the explicit expression of C r Eq. (|5.6p in the master 
formula Eq. f|3.33j) : 



a QCD (r, M 2 ) = a SCET (r,M 2 ) + a 2 s (M 2 ) + TW + ^'k) 



5=0 
(5.8) 



where 



J T Z % 

1 1 

n + £ n! 



oo 

(1 - rfA(0 £ — --(1 - r)VSa T (r, M 2 ); (5.9) 



n=0 
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4cet(t, M 2 ) = ^-a SCET (r, M 2 ), (5.10) 

up to corrections suppressed by powers of 1 — r, as shown in Appendix [Bl 

Equation (|5.8j) provides the sought-for explicit comparison of the QCD and SCET 
results at the level of hadronic cross-sections. Note that the non-convolutive nature of the 
SCET result implies that the generating function for the correction term is now given by 
the function S(r, £), which depends on the parton luminosity, rather than by the universal 
function K(z,£) Eq. (|577|) . 

In order to understand the correction term in Eq. (|5.8p . we note that, with the choice 
of /i s Eq. (|5.1|) . we get 



1 1 



c^(t,0 = e-^A(C) V -^(1 - r)V&Lr(r, M 2 ), (5.11 

^— » rj. 4- > n.i 



n=0 



n + ^ n\ 



so the dependence on (1 — r)^ cancels. It follows that £ derivatives acting on c^S(r, £) do 
not induce any extra logarithmic enhancement, other than that of <t scet (t, M 2 ) itself: 



° a 



n=0 

where the constants C n are r-independent: 



^qcd(t, M 2 ) = a SCET (r, M 2 ) + a 2 s (M 2 ) £ ^(1 - r)V^ ET (r, M 2 ), (5.12) 

n! 



3 ^ u d£ 3 8 ^ 2 z d£y n + e 



(5.13) 



2. . „ 7T 2 



Co = -gCa^iA) ~ ^A 2 , (5.14) 

A^o/vr 2 2\ A 2 2D 2 
C n = 2l~T~2 + ' n>0. 5.15 

Therefore, up to order a 2 , the correction term is just a 2 (M 2 ) times a linear combination 
of derivatives of a SCET with respect to ln(l — r). It follows that the correction term is at 
most of order 

a|xa^ln 2fc (l-r) xln p (l-r) = a^ln 2/l+p - 4 (l-r); h = k + 2, (5.16) 

where terms of order agln 2fc (l — r) are due to the coefficient functions, while terms of 
order ln p (l — r) are due to the parton luminosity. 

In other words, at order a", terms ln fc (l — r) in the SCET and QCD result coincide 
if2n — 3 + p < k < In. There are now various possibilities. If we simply neglect all 
logarithmic enhancements from the parton luminosity, i.e. if we set p = 0, then we conclude 
that the SCET and QCD results differ by terms which are NNLL according to the QCD 
counting Tab. [H but N 3 LL correction according to the SCET counting Tab. [2j Hence we 
conclude that, neglecting logarithmic enhancements from the luminosity, the SCET result 
does reproduce the QCD result to NNLL accuracy, albeit with the less accurate SCET 
definition of what is meant by NNLL. However, if a logarithmic enhancement from the 
luminosity is present, this is no longer the case, and the discrepancy can become arbitrarily 
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large (i.e. even at the leading log level) by just increasing the value of p. Note that this 
in particular means that with this choice of soft scale it is not possible to upgrade the 
accuracy of the SCET expression, as given in Tab. [21 to that of the QCD expression, as 
given in Tab. [TJ because at each logarithmic order the SCET expression differs from the 
perturbative QCD results by terms which, though consistent with the accuracy of Tab. [2 
spoil the higher logarithmic accuracy of Tab. [TJ 

One may then ask whether logarithmic enhancements due to the PDF are expected 
to be present, and whether they should be counted. Because the P qq and P gg splitting 
functions behave as P ~ (i- x ) + as 00 ~ * ^' con t r ibutions to all parton distributions f(x) 
which are enhanced by ln(l — x) terms will always be induced by perturbative evolution. 
In a parton distribution evaluated at some reference scale Qq these terms will be accom- 
panied by powers of a s (Qg), which is not small if the reference scale is taken as some low 
"initial" scale. Hence, in general, one does expect logarithmically enhanced contributions 
to PDFs, unless one wishes to make some fine-tuned assumption about the PDF itself, 
which can only hold at one single scale. The second question is whether these terms should 
be included or not in the power counting Eq. (|5,16p . This is a question which cannot be 
answered on the basis of first principles. Two relevant observations here are the following. 
First, once one substitutes any explicit expression of the parton luminosity in the expres- 
sion Eq. (|5.8p for the difference between the SCET and QCD result, there is no way to 
separate what comes from the luminosity and what comes from the coefficient function, 
because the SCET expression is not factorizable. Hence, in order to discard the luminosity 
logs from the power counting ones has to invoke the explicit SCET expression Eq. (|2.12p : 
in other words, one must argue that the SCET expression contains more information than 
that which is contained in the order-by-order perturbative result. The second observation 
is that in practice these correction terms may be parametrically large in realistic situ- 
ations, and they may lead to significant discrepancies between the predictions obtained 
using <7 SCET (r,M 2 ) or ct qcd (t,M 2 ). 

5.2 All orders 

The fixed O(o£) computation of Section [5.11 can be easily generalized to all orders. First 
of all, we note that the argument is based on the observation that to 0(a^) the correction 
term Eq. (|5.8p can be expressed as a series of derivatives of the function c^S(r, £) with 
respect to £ , but these do not lead to an extra logarithmic enhancement beyond that which 
is already present in S(r, £). The argument of Section [5. II would thus hold, to NNLL but 
to all orders in a s , provided only the correction term in Eq. (|5.8p was a series of series of 
derivatives of the function c^S(t, £) with respect to £ to all orders in a s . This is true if 
and only if S r depends on fi s only through powers of In with c given by Eq. (|5.5p . 
Now, we observe that the generic term in <S r , Eq. (|3.32|) . has the form 

4 7r«s(^)(m^-ln|^) =jf f<(^)(ln^-ln*J ,(5.17) 

with n > 2 and m = 0, 1, 2. This is not a function of In -ft only, because of the dependence 
of a s on jA s . In order to generalize the argument to all orders we must thus study this 
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dependence. Note that because 




(5.18) 




) 



rn 



m + 1 



1 




(5.19) 



For n = 2, this term contributes to Eq. (|5.8p an order-a| correction. This gives a series 
of extra contributions to the correction term, on top of those whose order was given in 
Eq. (|5.16p . which are at most of order 



a|ln(l-r) xa^ln 2fc (l-r) xln p (l-r) = a s l ln 2/l - 5+p (l-r); h = k + 3 (5.20) 



while higher-order terms are even more suppressed. 

The power counting which ensues from Eq. (|5.2(jp is the same as that of Section 15.11 
neglecting logarithmic enhancements from the luminosity (i.e. if p = 0), the SCET result 
does reproduce the QCD result to NNLL accuracy, but with the less accurate SCET 
definition of what is meant by NNLL. If a logarithmic enhancement from the luminosity 
is included, the QCD and SCET results differ, with the discrepancy appearing at any 
desired logarithmic order (including leading log) if the enhancement of the luminosity is 
sufficiently strong. 

It is interesting to observe that a different power counting might be appopriate in the 
phenomenologically relevant case in which the hadronic r is actually far from threshold, 
yet threshold resummation effects are non-negligible because the partonic center-of-mass 
energy is lower than the hadronic one, as discussed in the beginning of Section [3l In this 
case, it might be appropriate to simply take fj, s as some numerical constant, and compare 
directly the SCET and QCD expressions of the coefficient function C SCET and C QCD through 
the master formula Eq. (|3.33|) . But if r is far from threshold, then fi s Eq. (|5.1|) is of the 
same order as the hard scale M, As a consequence, in this case C r (N, M 2 , /i 2 ) Eq. ()3.3ip 
manifestly starts at NLL order, as confirmed by inspection of Eq. ()5.4|) . which gives 



One must therefore conclude that this class of NLL terms are resummed, through the 
leading-log function g\ in Eq. (|2.9|) . by the QCD result C QCr> (z, M 2 ) but are not resummed 
at all in C SCET (z, M 2 , fi 2 ). This problem may be alleviated through generalizations of the 
choice Eq. (|5.ip such as those proposed in Ref. [25] (and mentioned above at the beginning 
of Section [5]), whereby one rescales fj, s Eq. (|5.ip by a factor (determined for instance using 
scale-optimization methods), because they lead to smaller values of /j, s . These choices, 
however, do not affect the counting of logs, and it is therefore impossible to assess their 
impact in a comparison with standard QCD results, other than by numerical methods. 




(5.21) 
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6 Summary 



We have analyzed in detail the relation between the approach to threshold resummation 
based on perturbative factorization of Refs. [2H6|, and the SCET approach of Refs. [23H26] . 
with the main goal of exploring the viability, both theoretical and phenomenological, of 
the SCET prescription to treat the divergent nature of the perturbative QCD expansion 
in the soft limit. By deriving a master formula which connects resummed results in these 
two approaches, we have shown that the way they are related depends on the choice of 
soft scale in the SCET expression. We have explicitly performed calculations up to next- 
to-next-to-leading logarithmic accuracy, though it is easy to convince oneself that the 
structure of our master formula holds to any logarithmic order. 

We have shown that if SCET resummation is performed in Mellin space, then it coin- 
cides with the standard perturbative result. The SCET and QCD results then have the 
same accuracy, and are both beset by the problem of the divergence of the perturbative 
expansion. With this (partonic) choice of soft scale SCET and QCD provide alternative 
ways of deriving the same resummed result. 

If SCET resummation is performed in momentum space, as advocated in Ref. |23j, 
the SCET and QCD results differ by a non-universal term, which depends on the parton 
luminosity (explicitly given up to 0(a 2 s ) in Eq. (|5.4p ): the SCET approach separates 
off the series of divergent contributions which is then contained in this term, with the 
SCET resummed result now given by a convergent perturbative expansion. The price to 
pay for this is fourfold. First, because the difference term is non-universal, it may spoil 
the logarithmic accuracy of the resummed result depending on the parton luminosity. In 
particular, if the parton luminosity contains logarithmically enhanced contributions (as it 
generally will, based on its behaviour upon QCD evolution), the difference term may enter 
at any logarithmic accuracy (including at the leading-log level), unless one decides that 
logarithms coming from the luminosity should not be included in the power counting. Note 
however that, because this correction term is not factorized, there is no way of actually 
isolating the logs that come from the luminosity, other than to assume that the luminosity 
does not contain any. 

If this problem of non-universality is neglected, the SCET and QCD results are equiv- 
alent, however only by redefining the logarithmic accuracy to be always by one power 
lower, according to the counting of Tab. [2] rather than the more accurate perturbative 
QCD counting Tab. [TJ Hence the second price to pay is that the logarithmic accuracy 
of the SCET result in this case is always lower by one power of log, to all orders in a s . 
Third, while perturbative QCD resummation prescriptions such as the minimal |1L)| or 
Borel [12,13J prescription introduce corrections to the perturbative result which are power 
suppressed or more, the SCET prescription introduces a deviation which is only logarith- 
mically suppressed. And finally, the power counting and suppression in the SCET result 
must be done at the level of the hadronic scale 1 — r, while in QCD it is done at the level 
of the partonic scale 1 — z. In many cases of physical interest [27] it may turn out that 
the latter is small even when the former isn't: in these cases the QCD counting will be 
more accurate. 

It will be interesting to investigate the phenomenological implications of this state of 
affairs. Our result enables such an investigation, by providing a closed-form expression for 
the difference between the SCET and QCD results. 
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A Mellin transforms 



We collect here some useful results on Mellin transforms, while referring to Section 2 of 
Ref. [6] and the appendices of Refs. [13j[27] for a fuller treatment. 

The Mellin transform which are necessary for the computation of resummed terms, 
such as the exponent of Eq. (|2.5p , can be performed using 



dz z 



N-l 



ln p (l - z) 
1-z 



E 

k=0 



r (*)(!) d k f 1 - 1 /" \n p (l-z) 



k\ dL k 



dz ■ 



i-, +0 UI < A1 > 



1 



where L = In Eq. (127TO . 
At NNLL 



dz z 



N-l 



F (ln(l - z)) 
1-z 



d 1 



2 



7T 2 \ 



2 1 



6 J dL 2 



l - l ' N dz F(\u(l-z)) 



1 - z 



r N :i LL -i- O ( i 



(A.2) 



(where 7 = — r'(l) is the Euler constant) for any function F(£) which admits a Taylor 
expansion around I = 0. Now, 



d 7 2 <i 2 



1 -73T- + 



dL 2 ciL 2 



L p = L p - 7 pL p - L + 



-1 , p(p-l).,2 LP -2 



= (L- 7 ) p + 0(L p - 3 ). 
Hence, to NNLL accuracy, Eq. (jA.ip can be written in the equivalent form 



(A.3) 



dz z 



N-l 



F(ln(l - z)) 
1-z 



1-1/iV F(l n (i_z)) ^2 d 2 /-i-i/JV F(ln(l-z)) 



dz ■ 



1 - 2 



12 dL 2 



dz ■ 



1 - z 



(A.4) 



where 



A = Ae 7 ; ln^ = L- 7 . (A.5) 
N 

An essential ingredient in the discussion of Section [5] is the inverse Mellin transform of 

(A.6) 

where c is a constant. In order to compute it, we start from the identity 

1 d r ' 



ln r ' 



A d£ 



-A(0 t dzz N - l \^- 1 - 



(A.7) 



where A(£) = pT^y- Eq. (|A.7|) should be (and usually is) written in the form 



1 ^ n f 1 

i n n_L= —Am / dz^- 1 



ln«" 



n0) 



(A.8) 



5=0 
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so that the integral is well defined even when the derivatives and the limit £ — > are 
taken under the integral sign. This is not necessary for our present purposes. Rewriting 
Eq. (|A.7p with N replaced by N/c we obtain 



N ^'Jo z 

d£, n Jo z 



(A.9) 



after rescaling the integration variable z — > z c . The inverse Mellin transform of In™ -£ can 
now be immediately read off Eq. (|A.9p . 

B Convolutions 

In this Appendix we compute the integral 

£(r,£) = A(£) f'^af-)^- 1 -, (B.l) 

J rj- 2/ \ Z / Z 

where <r(r) = cr SCET (r, M 2 ) for simplicity, up to terms suppressed by powers of 1 — r. 
Using 



lni = l-z + 0((l-z) 2 ) (B.2) 



we find 



^)=a (5 )/;^(i-^-v(£ 

n=0 ' J T 

Expanding 1 / z n+l in powers of 1 — z we see that the integral is a sum of terms proportional 
to 

(1 - rf +m ; m>n. (B.4) 
Hence, the derivatives a^(r) appear in S(t, £) multiplied by (1 — r) m , with m > n. Now 

_ d(J ( T ) 



(l-r)a w (r) 



dln(l-T) 



(l_ r )V 2 )(r)- ^ (T) 1 d2(7{T) 



dln(l-r) dln 2 (l-r) 

(B.5) 

which means that (1— r) m cj ( - n - > (t) with m > n is power-suppressed, and can be we neglected 
in Eq. (|B.3|) since we are only interested in logarithmically-enhanced contributions to 
£(t, £). Hence 

^( r i£) = — r)^ — — — h power-suppressed terms. (B.6) 
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